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The eigenvalue problem for one-dimensional Schrodinger equation with the rational potential is 
numerically solved by the operator method. We show that the operator method, applied for solving 
the Schrodinger equation with the nonpolynomial structure of the Hamiltonian, becomes more 
efficient if a nonunitary transformation of the Hamiltonian is used. We demonstrate on numerous 
- - - examples that this method can handle both perturbative and nonperturbative regimes with very 

■ high accuracy and moderate computational cost. 
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O ■ I- INTRODUCTION 



PACS numbers: 03.65.G,03.65.F 



In this paper we generalize the operator method for solving the Schrodinger equation (SE) with a non- 
polynomial structure of the Hamiltonian. As an example of such a problem we consider one-dimensional 
Hamiltonian with the rational potential 

■ 

' The problem chosen is of interest in different areas of physics, i.e. laser physics! and non-linear quantum 

CN| \ field theory^. It can be also exploited for verifying new nonperturbative methods, which allow us to find 

■ approximate solution of the Schrodinger equation in the whole range of Hamiltonian parameters. Such 

methods can be used for treating many physical problems of current interest. The test problem chosen 
is very suitable for these purposes because the SE 



^ + Vl{x)~E 



^j) = 0, (2) 



' can be solved exactly if some algebraic relations between parameters g and A are assumedSiii^. The exact 

^ I solutions allow to find out the accuracy of approximate methods. 

O^. The Schrodinger equation with the potentials Vi{x) and V2{x) has been treated by different approaches 

based on finite-difference schemes, variational and asymptotic methodsSi&^iiSiiiiiSiiSi. To the best of our 
- ^ . knowledge there is only one worki^, where the case L > 2 was studied in the wide range of the parameter 

' a = g/Xhy very efficient method of finite-differences. The case L = +3 has been studied for g <C 1, A <C 1 

^ I by using the Pade-amjroximant and hypervirial methods!^. Three-dimensional generalization of Eq. (01 

can be found in Ref. [l^ 

In the present paper the SE Eq.lj^J is solved by the operator method (OM), developed by Komarov 
and Feranchuki^. The OM has been successfully applied to different problemsi2iiSii2iS£, but the most 
problems considered had a polynomial structure of the Hamiltonian. The polynomial structure allows 
one to calculate matrix elements of the SE equation analytically in the Fock basis^^. But this simple 
analytical scheme can not be directly applied to the nonpolynomial equation Eq. Q . We can still perform 
the OM calculation&l^, but some extra numerical scheme must be exploited to calculate matrix elements. 
This problem can be resolved by the use of the nonunitary transformation C 

1^-)=^!^), C={l+gx^), (3) 

which leads to the following expression for Eq.® 



y + y + .9(^ + y)+Ax^^-(l + .9x2)£; 



l^> = 0. (4) 



So we transformed the nonpolynomial operator equation Eq. Q to the polynomial one Eq. I^J , where the 
operator at the left side is nonhermitian with respect to the scalar product 

(cj^l^P) = / dx^ixm^). (5) 
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The idea of using nonunitary transformations was recently suggested in Ref. 21. It was also developed in 
Ref. [23 as nonhermitian technique of canonical transformations. The basic idea is that the commutation 
relation [x,p] = i can be preserved by any unitary transformation, but the same is also true if we make 
use of the similarity transformation X — CxC^^, P — CpC~^. Below we show that the polynomial 
operator equation Eq.Q) is very easy to handle within the framework of the operator method in order 
to find both the exact numerical and approximate analytical solutions. The OM allows us to treat the 
problem in the wide range of Hamiltonian parameters g and A, where the SE Eq.Q has a discrete energy 
spectrum, i.e. when -00 < L < +1, |A| < +00, g > 0; L > +2, A > 0, g > 0; L = +2, A > -g/2, g > 0. 
Below we present results of our nonperturbative calcutations for four cases I < L < 4. 



II. ITERATIVE SCHEME OF THE OPERATOR METHOD 



In this section we apply the operator method for solving the equation Eq. Q using the scheme developed 
in Ref. 16. If we introduce the annihilation operator a(w), the creation operator a+(ti>) and the excitation 
number operator n(cj) — a'^{uj)a{uj) 

^'^i^) ^ ^=i^^ - ^P)^ a(w) = — i=(tja; + ip), (6) 



which satisfy the commutation relations 

[a{Lu), a'^ [uj)] ~ 1, [a{LL>),n{uj)] = a{uj), [a'^ (uj),n{Lu)] = —a^ (uj), (7) 
one have for Eq.l^J 

(£' - EC"M=0, (8) 
£' = (i- + ^)(2.+ l) + |j(2.^ + 2.+ l) + |(2.^ + 2.-l) + |(-l-l) 



X (a+4 + a*) + 



/I ^ 3.9 , .9 \ / +2 , 2n 9+2 



+ |(a+2-a2) + A.T2^ (9) 

C" = l + ^{a+^+a^ + 2n + l), (10) 
2uj 

where x^^ is written below for the four particular cases 1 < L < 4 

x' = l-(2n + l) + i-(a+2 + a2), (11) 

x^ = A(i + 2n + 2n2) + JL(a+4 + a4 + (4n-2)a+2 + (4n + 6)a2), (12) 
4(jj^ 4(jj^ 

— ^(3 + 8n + 6n^ + 4n^) + ^(0+*^ + + (6n - 9)a+'' + (6n + 15)a'^ 
+ 15(1 - n + n^)a+^ + 15(3 + 3n + n^)a^), (13) 
^ (105 + 28071 + 35071^ + 140^3 + 7071*) 



x^ 



16^4 



^ -(a+« + a* + (871 - 20)a+^ + (871 + 28)a^ 



16u;4 

+ (98 - 8471 + 287i2)a+'' + (210 + 14071 + 2871^)0" 

+ (-84 + 19671 - 8471^ + 5671^)0+2 + (420 + 5327i + 25271^ + 5612^)0^). (14) 

Any degree of x'^ can be derived by the use of a recursive relation between a;^^ and x^^^^. The normalized 
eigenvectors \n,uj) of the excitation number operator 71 are 

\n,u:)=^-^-^\Q,Lo), a{uj)\0,u;) = 0. (15) 
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In contrast to the equation Eq. (01 we can calculate the matrix elements of the equation Eq. ((Sj) exactly 
in the Fock basis Ea. H15|) . The eigenvector \ip) can be represented in the Fock basis Ea. ljlSII as 

|^)=^Cpb,c.), (16) 

with the coefficient Cp to be defined by the equation Eq.Q. The arguments of a'^(uj), a(w), n(w) and 
\n,uj) will be omitted in order to simplify all expressions. Substitute Ea. l|lt)|) in Eq.® and find its 
projection on bra-vector {k,uj\ 

Ck{k\C'\k) - Ck{k\C"\k)E + Y,Cp{k\C'\p) ~ EY,Cp{k\C"\p) = 0. (17) 

p^k p=jtk 

The system of the linear equations Ea. (|17|l can be solved by the iterative method^Si 



{k\a\k) - Ei'-^\k\L"\k) 

{n\C'^n) ^YrpZ\.^p^rA'-'\n\a'\py ' ' 

(s) (s) 

where is the coefficient of the eigenvector expansion and En is the energy of n-level calculated in 
s-th iteration; m = 2L; s = 0, 1, . . . , Smax- The results of numerical calculations of the energy eigenvalues 
for the ground (n = 0) and first excited {n = 1) states in the wide range of values of g, X are summarized 
in Tables 1-4. For convenience of the comparison with the results in Refs. Il2ll4l we present our results 
in the form 2A, 2E. As is clear from Table 1 and Table 2, where some exact eigenvalues are given, the 
accuracy of the operator method is very high lEn"""^"^ — £'„]/£'„ < 10^^^. But the case g ^ 1, A ^ 1 
needs the increase of the iteration number Smax- So it is important to make use of the parameter u, 
which allows us to speed up the convergence of the iterative scheme. Since the exact eigenvalues of the 
Hamiltonian do not depend on the choice of uj, the following condition has to be satisfied for exact energy 
eigenvaluesi^ 

0. (19) 



duj 



Because En is equal to E^\uj), a good accuracy can be achieved at the extremum point ujextr of the 
function (w) (see Fig. 1). Notice that this statement is only valid in the case of large iteration 

number. As soon as the iteration number is relatively small En""^'^\Lu) is described by an oscillating 
function crossing E = i?„ (exact) for some values of uj (see Fig. 1 and Ref. ITgl) . In the case of large 
iteration number the magnitude of oscillations is very small and goes to zero quickly (see Fig. 1). 

We emphasize that even for small number of iterations the points of extremum are near the exact 
energy eigenvalue and belong to the range of stable convergence of the iterative scheme. They can be 
chosen as starting values for lo. As is clear from Table 2 the values of Smax and u must be increased for 
excited states. Large values of g and A require the increase of u as well. Such behavior is typical for all 
potentials Vl{x). Different criteria for choosing u can be also found in Refs. T^Tz. 19. 

We'd like to notice that the normalized eigenfunctions '^n{x) of the original SE Eq.© can be easily 
calculated 

^n{x) = « 1(1 +gx^)5c^-)(a:|fc,c.) (20) 



k=0 



N' = (V'nIV'n) « / dx{l+gx') 



Lfc=o 



(21) 



The polynomial structure of the operator equation Eq.Q) also allows us to obtain some analytical ex- 
pressions for the energy eigenvalues with the explicit dependence on all parameters. Such expressions 
can be found by making a limited number of iterations with any system of computer algebra like Maple 
or Mathematica. 
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III. CONCLUSION 



The Schrodinger equation with the rational potential was solved by the operator method. The energy 
of the ground state and the first excited state were calculated with very high accuracy without any 
limitations on the parameters of the Hamiltonian. The numerical solutions are in excellent agreement 
with the exact solutions known for some values of g and A. The criteria for choosing the convergence 
parameter uo was analyzed. So we have shown that the spectrum of the original nonpolynomial Schrodinger 
equation Eq. is identical to the one of the polynomial operator equation Eq. I^J . The results obtained 
demonstrate that the operator method is very efficiently for solving the problems with the nonpolynomial 
structure of the Hamiltonian if the nonunitary transformation is used. The application of the operator 
method to the three-dimensional nonpolynomial Hamiltonian is straightforward. 
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TABLE I: Energy eigenvalues for Vi{x). 



2\/g 100/10 10/10 10/100 

2Eo{OM) 5.79394230019270 1.58002232739150 1.08406333549441 
2Ei{OM) 11.5721967757092 3.87903683088257 3.09831699508894 
2X/g -2.494002/0.499 -1.585842/0.339 -0.0402/0.01 

2^0 (OM) 0.002000000000000 0.322000000000000 0.98000000000000 
2£;o (exact) 0.002000000000000 0.322000000000000 0.98000000000000 



TABLE II: Energy eigenvalues for V2{x). 



2X/g 


0.1/1 


0.1/0.1 


1/0.1 


2£;o(OM) 


1.02514716380525 


1.05529770725788 


1.36059173241772 


^/ Smax 


3.0/69 


1.74/20 


2.90/22 


2£;i(0M) 


3.09577659167734 


3.24865072344421 


4.49215757165382 


W/Smax 


2.9/70 


1.56/20 


2.82/22 


2X/g 


100/10 


10/10 


10/100 


2Eo (OM) 


2.90257776978742 


1.35977466215783 


1.04797115328336 


^/ ^rnax 


16.6/115 


11.0/200 


27.3/1400 


2Ei (OM) 


9.21956081849447 


4.15798607191020 


3.14544411596743 




16.3/125 


10.4/210 


25.7/1700 


A/g 


0.001/0.02 


0.005/0.1 


0.01/0.01 


Eo{OM) 


0.500712692243485 


0.503010339434091 


0.507093239851910 




0.5007126922434854 


0.503010 


0.507093239851 


Ei{OM) 


1.503492319497928 


1.51392272439619 


1.53457040870676 


Ei^) 


1.5034923194979289 


1.51392 


1.534570408 


a/g 


0.5215/0.2607543752089 






2^0 (OM) 


1.18657366612601 






2£;o(exact) 1.18657366612601 







TABLE III: Energy eigenvalues for V3{x). 



2X/g 


100/10 


10/10 


10/100 


2Eo (OM) 2.32812938769652 


1.36685074612832 


1.06484795071250 


2Ei (OM) 8.25858912303721 


4.57644754780293 


3.30561773097224 


A/<? 


io~Vio-3 


lO'^'/O.Ol 


5 X 10"''/2.5 X 10"^ 


Eo{OM) 


0.500186311611817 


0.500180733817835 


0.500916657050025 


EoC') 


0.5001863116118168 


0.5001807336178352 


0.5009166570 


2X/g 


10/0 


1000/0 


10000/0 


2i;o(OM) 


2.20572326959563 


6.49235013232967 


11.4787980422645 


2Eol^) 


2.205723269595632 


6.492350132329672 


11.47879804226454 
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TABLE IV: Energy eigenvalues for V4{x). 



2\/g 100/10 10/10 10/100 

2Eo (OM) 2.14407075597610 1.41752810562730 1.10864509823789 

2Ei (OM) 3.59432196739250 4.96949962368201 3.59432196739250 

2X/g 10/0 1000/0 10000/0 

2£'o(OM) 2.11454462194213 4.94948744003274 7.77827221431110 

2Eo^) 2.114544621942129 4.949487440032743 7.778272214311099 




1.0-1 1 1 1 1 1 U! 

1 2 3 4 5 

FIG. 1: The ground state energy 2E^''^ as a function of for L = +2, 2£;o(exact) = 1.18657366612601, 2A = 
0.500008389662453, g = 0.260754375208969. 



